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Abstract 

A generalization of the Polyakov-Koba-Nielsen-Olesen scaling law of the mul- 
tiplicity distributions P{n, s) is developed. It states that a suitable change in 
the normalization point of P{n, s) compensated by a rescaling can restore data 
collapsing onto a universal curve if the original scaling rule is violated. We 
show that the iteratively executed transformation of P{n, s) can be viewed 
as varying the collision energy. The e~^e~ and pp multiplicity data at top 

energies are found to exhibit a fixed point property of the iteration. 
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This year is the 25th anniversary year of the influential work of Koba, Nielsen and 
Olesen concerning the asymptotic scaling behavior of multiplicity distributions |1|]. They 
put forward the hypothesis that at "sufficiently high" energies s the probability distributions 
P{n, s) of detecting n final-state hadrons in a certain collision process should exhibit the 
scaling relation 



where {n{s)) is the average multiplicity at collision energy s. According to Eq. (|l|), P(n, s) 
is a homogeneous function of degree —1 of n and {n{s)). The homogeneity rule states that 
the multiplicity distributions are simple rescaled copies of the universal function ip, i.e., 
the change of collision energy s amounts only to a change of scale in the shape of P{n,s). 
The above prediction caused immediately a great deal of activity in the experimental and 
theoretical analysis of multiplicity distributions which remained immense during the past 
25 years. It should be mentioned here that Polyakov arrived at the scaling law Eq. (|1]) 
already in 1970 postulating a similarity hypothesis for strong interactions in e+e^ —>■ hadrons 
annihilation 0. Despite of this fact, the phenomenon became known as KNO scaling in the 
high-energy physics community. 

The main goal of the present Letter is to extend the domain of data collapsing behavior of 
P{n, s) beyond the scaling relation Eq. (|l]). In particular, we propose an iterative procedure 
capable of finding homogeneity rules for P{n, s) if the scaling hypothesis Eq. (|l]) is violated. 
We shall consider in detail the fixed points of the iteration, their domains of attraction and 
the constraints they put on the asymptotic scaling form of P{n, s). Our approach is inspired 
by the renormalization group methods in the theory of critical phenomena. In the beginning 
of the paper we will review some basic results of Ref. 0] which should be considered as the 
first part of the present work; it is cited as Paper 1 in the foUowings. 

We start with some elementary properties of the scaling function ip{z), z denoting the 
scaled multiplicity n/ (n(s)). Finite energy discreteness effects at moderate (n(s)) are usually 
taken into account via Poisson transform which is discussed in Paper 1. The continuous 
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probability density ip{z) fulfills the normalization conditions il){z)dz = zip{z)dz = 1. 
The latter one, yielding the constraint (z) = 1, defines a second properly normalized scaling 
function: ip'lz) = ziIj{z) = nP{z). Obviously, the moments of ip{z) coincide with the 
normalized moments of P(n, s), 



oo 



C, = i z'^mdz = (2) 







They are independent of collision energy s if the scaling hypothesis Eq. (|l|) holds valid. The 
moments of ^ip'^z) are given by C'^ = /o°° z'^ ip\z)dz = C^+i i.e. the difference between the two 
scaling functions amounts to a shift in the ranks of their moments. 

Let us now recapitulate a scaling argument first presented in Paper 1. Assume that 
violation of the original scaling hypothesis Eq. (|l]) is observed and we measure energy de- 
pendent "scaling" functions ip{z,s) and ip'{z,s)- A key point to the later developments is 
the observation that ip'{z, s) enables one to perform a rescaling of type Eq. (|l|) yet again to 
arrive at data collapsing behavior. The modified scaling hypothesis reads 

where the subscript of the new scaling function refers to the first step of an iterative sequence 
of transformations that can be performed on P{n,s). For quantities corresponding to the 
0th step the subscript will be omitted. Our new scaling variable is zi = z/C[{s) = z/C2{s) 
and the first-iterate moments Cg^i = z^ ipi{zi)dzi are related to the original ones through 

Expressing ipi{zi) in a similar manner we get ipiizi) = C2{s) n P[z/C2{s)] which can be 
rewritten according to 

It is seen that obeys the same structure as the original scaling function ■^/'(z), namely, 

iIji{zi) = {n{s))i ■ Pi[n/ {n{s))i]. The first-iterate multiplicity distribution Pi(n, s) and its 



moments are obtained by changing the normahzation point of P{n,s) from the 

0th moment to {n{s)) and rescahng by {n{s)) to maintain the overall normalization. The 
distribution Pi (n, s) is known in the mathematical literature as the first order moment 
distribution of P{n, s). 

If ipii^i) is not independent of collision energy s, the transformation rule described above 
can be repeated iteratively until the appearance of data collapsing onto a universal scaling 
curve. The details are presented in Paper 1, here we recall only the final result. In the ith 
step of the iteration the connection between Pi{n, s) and the original multiplicity distribution 
P{n, s) is provided by 

P^(n,s) = —^P{n,s) (6) 

i.e. the normalization point of P{n,s) is changed to (n*(s)) and a rescaling by {n^{s)) is 
made to preserve proper normalization. Thus the necessary condition of performing the ith 
iteration step is the existence of the moments of P{n, s) up to ith order. Pi{n, s) is the mo- 
ment distribution of order i of P{n, s). Data collapsing of Pi{n, s) onto the scaling function 
= ■ occurs whcu the normalized moments at the previous itera- 

tion step exhibit the "monofractal" type behavior Cq,i_i(s) oc [C2,j_i(s)]'^^^ with constants 
of proportionality independent of collision energy s, see Eq. (^). It is worth mentioning 
here that for discrete probability laws such as P{n, s) the factorial moment distributions, 
involving factorial powers of n, arise more naturally as discussed in Paper 1. 

It is also argued in Paper 1 that the iterative procedure acting on P{n, s) bears some 
similarity with the renormalization group (RG) transformations. They correspond in many 
important applications to a change in the norm of the parameters characterizing a physical 
system . On the one hand, Eq. (|^) can be viewed as a Gell-Mann - Low type relationship 
in which a multiplicative transformation (here multiplication by n*) is compensated by a 
rescaling and a suitable change in the renormalized parameters (here in the moments). Since 
the operation Pi{n, s) Pi_i(n, s) does not exist the iteration steps constitute a semigroup. 
In statistical physics the Kadanoff - Wilson type RG transformations perform a systematic 



reduction in the number of degrees of freedom via e.g. spin decimation. This makes possible 
to ehminate the small-scale fluctuations from the problem which are irrelevant to critical 
point behavior such as the homogeneity of thermodynamic functions [Q. In the iterative 
procedure described above the elimination of small-scale fluctuations corresponds to the 
elimination of the low-order moments of P{n, s) via moment-shifting and rescaling until one 
observes a homogeneity rule of type Eq. (|l|) for Pi{n, s). 

The suggested analogy between the transformation Eq. (|) of multiplicity distributions 
and the RG methods in field-theory and statistical physics may seem too remote at first 
glance. But there is a very close relationship between RG ideas and certain concepts of 
probability theory ||5|Q which naturally fits into our approach to asymptotic multiplicity 
scaling. Let us demonstrate it by considering some properties of Eq. (|]) in more detail. 

Property #1: Form-invariance under size-biasing 

If we view the Gell-Mann - Low or Kadanoff - Wilson version of the RG as a trans- 
formation acting on a probability distribution, it is a transformation that does not change 
the form of the distribution 0. Weighting a probability law according to Eq. is known 
in statistics as size-biasing of order i. Members of the log-exponential family of distribu- 
tions, such as the beta, gamma, Pearson type V, Pareto and log-normal to mention but a 
few, are form-invariant under size-biasing, i.e. they retain functional form and only their 
parameters are affected 0. Therefore if P{n,s) belongs to the log-exponential family, our 
iterative procedure is in accordance with this particular aspect of the Gell-Mann - Low and 
Kadanoff - Wilson type RG transformations. 

Property #2: Fixed points and automodel distributions 

A probability distribution is called a scaling- or automodel distribution if it is invariant 
under the action of the RG, or, in other words, if it is a fixed point of the RG transfor- 
mation [^]. In our case automodel distributions are those members of the log-exponential 
family which are affected by size-biasing through a scale-change. A theorem of Ref. states 
that size-biasing amounts to a change of scale in the original distribution if its normalized 
moments have the form Cq = C^^'' This is a well-known property of the log- normal 
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law whose parameters are transformed by Eq. according to u ^ u + ia"^ with u and cr^ 
being the mean and variance of the unbiased distribution. A log- normal s) displays 
fixed point behavior in the following manner: although the transformation Eq. (|^) changes 
its first moment, this change is scaled out by constructing ipi{zi) and one arrives at an un- 
affected scaling function at each iteration step. Therefore a log- normally shaped ipi^) is a 
fixed point of the iterative procedure. Let us emphasize that the log-normal law is not a 
unique fixed point of the iteration because the distribution is not uniquely determined by 
its moments. 

Property #3: Domains of attraction of the fixed points 

Besides finding the fixed points of the RG transformation, it is important to specify their 
domains of attraction. A domain of attraction is the set of initial probability laws which 
converge to a given automodel distribution under the action of the RG transformation 0. 
The domains of attraction are analogous to the universality classes of critical phenomena. 
Here we are interested in those probability laws which converge to the log-normal fixed point 
under the action of Eq. (^. Let us choose ipi^) to be the generalized gamma density 

V^(^) = ^A^V=-iexp(-[A.]^) (7) 

with shape parameter k > 0, scaling exponent yU > and scale parameter A restricted to 
A = r{k + l/fi)/r{k) by the normalization condition (z) = 1. In the limit fi ^ and 
/c — s> oo the scaling function given by Eq. (|^ converges to a log- normal ip{z) of variance 
0"^ in such a way that kfi"^ — > 0. The generalized gamma distribution is a member 

of the log-exponential family and thus it is form-invariant under size-biasing. The shape 
parameter k is changed by Eq. according to k —>■ k + i/fi whereas the scaling exponent n 
remains unchanged |[T^. The fi ^ 0, k —>■ oo log- normal limit can be achieved by requireing 
C2,i = C2 in the course of the iteration; the increase of k is compensated by a decrease of fi 
so that the second moment of iplz) is not affected by the variation of the parameters. Thus 
suitably standardized generalized gamma distributions constitute a domain of attraction of 
the log-normal fixed point of the iteration. 
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Property #4- inequality for the first moments 

The transformation rule Eq. is a special case of weighting a probability distribution 
f{x) with the non-negative weight-function w{x) according to w{x)f{x)/{w{x)). The first 
moment of the weighted distribution is greater or smaller than the first moment of f{x) 
depending on whether w{x) is a monotonously increasing or decreasing function of x. In the 
case w{x) = x, i.e. for size-biasing of order one, we have the inequality {x)i > (x) and in 
the general case > 

Consequences 

In many applications of renormalization group methods the transformations establish 
correspondences between physically different states of the same system. For example, in 
the Ising model the iterative repetition of the RG transformation can be viewed as varying 
the temperature. Similar correspondence can be established for asymptotic multiplicity 
scaling. According to the above inequality of the first moments, the iterative sequence of 
normalization point changing transformations acting on P{n, s) can be viewed as increasing 
the collison energy s. The asymptotic scaling relation Eq. (|T]), which states that the increase 
of s amounts to a change of scale in the shape of P{n, s), is therefore a fixed point property. 
An important consequence of our reasoning is the fact that the asymptotic scaling function 
can not be arbitrarily shaped: the above fixed point behavior is satisfied only by those 
probability laws which are equivalent in their moments to the log-normal distribution. 

In the light of the previous findings it is of interest to estimate the degree of deviation 
between the log- normal distribution and the shape of ip{z) at asymptotic energies. Making 
use of the generahzed gamma density Eq. (|^ the parameter pair {k, /i) is well suited to 
measure the departure from fixed point behavior corresponding to the fi ^ 0, k ^ oo limit. 
The shape of the asymptotic ip{z) can be reconstructed by fitting the Poisson transform 
of the generalized gamma distribution, let us call it HNBD for short, to the experimental 
data for P{n, s) available at the highest collision energies s. The HNBD was invented 
and developed in Refs. [|n],|T^ by the present author. The analytic form of P{n, s) can be 
expressed in terms of generahzed special functions: 
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where JC = k + n/ fi, = n\ T{k) 6^"^ and 6 = {n{s))rik) /T{k + of course the shape 

parameter k and scahng exponent fi may also depend on s. The functions H}'}(-) and iM/o(-) 
are particular cases of the Fox- and Wright hypergeometric functions, respectively |13| . The 
negative binomial distribution given by Eq. (P) is the /i = 1 marginal case of Eq. (^) for 
{n{s)) > k and of Eq. (p!OD for {n{s)) < k. The Poisson transformed log-normal limit (/i = 0) 
lacks a representation in terms of known functions. 

We have investigated two full phase-space data sets for P(n, s): the Delphi data at 
^/s = 91 GeV in e+e^ annihilations |T^ and the UA5 data at ^/s = 900 GeV in pp colli- 
sions [|T5[. For each data set we have performed several HNBD fits with different scaling 
exponents /x. The value of fi was varied in the interval 0.1 < ;U < 2 by stepsize Afi = 0.1. 
It is worth considering how the and the best-fit shape parameter k depend on fi. The 
trends are displayed in Fig. 1. The top right inset shows the variation of fit quality. As is 
seen the decreases monotonously towards fi = 0, approximately as an exponential, for 
both data sets. The shape parameter k increases according to a power-law as /i — > with 
slopes being the same for the two reactions (the estimated errors of k are too small to be 
seen). Let us recall that the convergence to a log-normal law of variance o"^ is such that 
/i — > and k ^ oo with kfi"^ l/o"^. Thus the /x-dependence of k allows us to estimate 
the value of cr by fitting k = [afi)"'^ to the data points. The fits are represented by the 
straight lines in Fig. 1. The quality of fits and the best-fit value of a are shown in the first 
row of Table I for each data set. The second row quotes the same numbers obtained in 
by fitting the Poisson transform of the log- normal distribution to P{n,s). All these results 
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suggest the same conclusion: the asymptotic scahng function V'(^); as can be guessed from 
pre-asymptotic multipUcity data, is log- normally shaped both in e^e~ annihilations and in 
pp collisions. Interestingly, this is just the fixed point behavior we would expect on the basis 
of the RG approach to asymptotic multiplicity scaling. 

Summarizing our results, we have developed an iterative procedure well suited to find ho- 
mogeneity rules of type Eq. (|l]) for the multiplicity distributions P(n, s). The quoted scaling 
law states that the energy dependence of P(n, s) is due entirely to the energy dependence 
of its first moment {n{s)). Thus rescaling P{n,s) by the average multiplicity according to 
Eq. (|l|) a universal scaling function emerges whose shape is independent of s. This scaling 
argument can be extended to the more general case when the s-dependence of P{n, s) can 
not be attributed to {n{s)) exclusively. To arrive at data collapsing onto a universal scaling 
curve, a modified scaling transformation is needed capable of eliminating the s-dependence of 
higher-order moments. An obvious solution is to "shift out" the moments of the multiplicity 
distributions up to the required order. This can be achieved by changing the normalization 
point of P{n,s) from the 0th moment to {n{s)), rescaling by {n{s)) to maintain the overall 
normalization and repeating the two-step transformation until the appearance of the scal- 
ing law Eq. (0) for the iterated distribution. In other words, we have the possibility that 
the increase of collision energy s amounts to a change of scale not in P{n,s), rather, in a 
moment distribution Pi{n,s) with i > 0. Accordingly, the phenomenon of data collapsing 
onto a universal scaling curve should be checked for i > as well in Eq. (^. During the 
past 25 years this was traditionally done only for i = 0. 

Besides a whole family of new scaling relations for the multiplicity distributions, Eq. 
provides a close analogy with certain RG ideas. In this respect there is a distinguished 
role of those probability laws which are form-invariant under the action of Eq. (|^) and the 
moment distributions Pi{n, s) are simple rescaled copies of P{n, s). These distributions are 
the only ones which may exhibit the asymptotic scaling behavior Eq. (|I]). It is known that 
all probability laws having the above property are equivalent in their moments to the log- 
normal law. Remarkably, the multiplicity data in e^e~ annihilations and in pp collisions 
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at the highest available energies indicate log-normality of the asymptotic scaling function 
tp{z) — just as one would expect on the basis of RG arguments. To decide whether this 
is purely accidental (as already happened with results obtained by related arguments |I6[ ) 
or is an intrinsic feature of multiparticle production, measurements at even higher energies 
are needed. The forthcoming pp data at ^/s = 1800 GeV at Tevatron are awaited with keen 
interest . 

This research was supported by the Hungarian Science Foundation under Grant No. 
OTKA-T024094/1997. 
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TABLES 



Data set 


a 


xVd.o.f. 


e+e", ^/s = 91 GeV 


0.199 ±0.001 


12.0/18 




0.201 ± 0.004 


32.4/24 


pp, ^/s = 900 GeV 


0.527 ± 0.003 


17.3/18 




0.538 ± 0.01 1 


32.7/52 



TABLE I. For each data set the first row displays tlic result of fitting the power-law k = (cr^)~^ 
to the /i-dependence of the HNBD shape parameter k. The fits are represented by the straight lines 
in Fig. 1. The second row quotes the outcome of Poisson transformed log-normal fits to P{n,s), 
see Ref. [12] for details. 
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FIGURES 




0.1 1 2 1^ 

FIG. 1. Variation of the best-fit shape parameter k of the HNBD as a function of scahng expo- 
nent iJ, for e+e" annihilations at = 91 GeV (open circles) and for pp collisions at ^/s = 900 GeV 
(solid squares). The straight lines are power-law fits, see the text and Table I for details. The inset 
in the top right displays the /x-dependence of corresponding to the HNBD fits. 
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